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Abstract 

We construct the fusion operators in the generahzed r*^^)-niodel using the fused L-operators, 
and verify the fusion relations with the truncation identity. The algebraic Bethe ansatz discus- 
sion is conducted on two special classes of r^^^ which include the superintegrable chiral Potts 
model. We then perform the parallel discussion on the XXZ spin chain at roots of unity, and 
demonstrate that the s^2-loop-algebra symmetry exists for the root-of-unity XXZ spin chain 
with a higher spin, where the evaluation parameters for the symmetry algebra are identified 
by the explicit Fabricius-McCoy current for the Bethe states. Parallels are also drawn to the 
comparison with the superintegrable chiral Potts model. 
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^In this paper, the term "XXZ spin chain" means a lattice model with the L-operator associated to the trigono- 
metric _R-matrix of the ice- type model. The "XXZ spin chain of higher spin" here is named as the "XXZ Heisenberg 
model with arbitrary spin" in [31]. 
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1 Introduction 



In this paper, we investigate the uniform structure about symmetry of various lattice models^, 
namely the generalized r^^^-model, the chiral Potts model (CPM), and XXZ spin chain of higher 
spin. Parallels are also drawn to differentiating the symmetry nature appeared in those models. The 
generalized r^^^ -model, also known as the Baxter-Bazhanov-Stroganov (BBS) model [3, 7, 9, 10, 24], 
is the six- vertex model with a particular field (see, [7] page 3), i.e. having the i?-matrix of the asym- 
metric six- vertex model, (see (2.3) in this paper). Note that the non-symmetric diagonal Boltzmann 
weights of this twisted i?-matrix differentiates it from the usual (gauge transformed) trigonomet- 
ric i?-matrix (see, for instance [11] (3.3) and references therein), and the usual construction of 
L-operator by employing the highest-weight representation theory of quantum groups or algebras 
fails in this context. However using the cyclic A^- vector representation of the finite Weyl algebra, 
one can construct a five-parameter family of L-opcrators for the A^-state generalized r'-^-'-model. 
When the parameters arc restricted on a family of high genus curves, called the rapidities, one 
obtains the r'^'^^-matrix of the solvable iV-state CPM. In [9] Bazhanov and Stroganov showed that 
the column-transfer-matrix of the L-operator in CPM possesses properties of the Baxter's Q-matrix 
(i.e., the chiral- Potts-transfer matrix) with symmetry operators similar to, but not exactly the same 
as, the SZ2-I00P algebra. Indeed the r^^^ -degenerate eigenvalue spectrum occurs only in superin- 
tegrable case [39], where the Onsager-algebra-symmetry operators in the quantum Hamiltonian 
chain [25] derived from the Baxter's Q-matrix provide the precise symmetry structure of the su- 
perintegrable r^^^ -model [40]. Due to the lack of difference property of rapidities, the study of 
CPM, such as the calculation of eigenvalues [4, 5, 34] and order parameter [8], relies on the method 
of functional relations among the t^J^- and chiral-Potts-transfcr matrices, consisting of (recursive 
and truncation) T-fusion relation, TQ- and QQ-relations. On the other hand, in the study of 
the "zero-field" six- vertex model ([6] Sect. 3) with the usual trigonometric i?-matrix, it has been 
extensively analyzed in [15, 18, 19, 41] that the degeneracy of the spin-i XXZ Hamiltonian occurs 
when the anisotropic parameter g is a root of unity with the extra sZ2-loop-algebra symmetry of the 
system [15]. Recently, the sZ2-loop-algebra symmetry appears again in the XXZ chain of spin-^^^^ 
at the Nth root of unity [36]^. Results about the root-of-unity symmetry of the six- vertex model 
were first discovered by algebraic Bethe ansatz and quantum group theory, later by the method 
of Baxter's Q-operator, and the degeneracy of eigenvalues was further found in the root-of-unity 
eight-vertex model [16, 17, 20]. By this effort, Fabricius and McCoy observed the similarity between 
the root-of-unity eight-vertex model and CPM, and proposed the conjectural functional relations 
of the eight-vertex model by encoding the root-of-unity symmetry property in a proper Q-operator 
of the theory, as an analogy to functional relations in CPM [21]. Along this line, the Q-operator 
incorporated with the s/2-loop-algebra symmetry of the spin-^ XXZ chain at the roots of unity was 
constructed in [42] where functional relations in the Fabricius-McCoy comparison were verified. By 

^AU the models discussed in this paper will always assume with the periodic condition. 

^The XXZ chain of spin--^^ at the iVth root of unity in this paper is phrased as the Nilpotent Bazhanov- Stroganov 
model in [36]. 
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this, we intend to make a detailed investigation about common features shown in those known mod- 
els, compare the symmetry structure, then further explore more unknown models in the scheme of 
functional relations. Among the function relations, the fusion relation with "truncation" property 
plays a vital role in connecting the T- and Baxter's Q-matrix. Usually one can derive the boundary 
fusion relation, i.e. the "truncation" identity, using the Baxter's TQ-relation for any Q-operator 
(not necessarily satisfying the QQ-relation) if a Q-matrix exists (see, e.g. [2, 29, 37]). Conversely, 
the validity of boundary fusion relation may as well strongly suggests the existence of Baxter's 
Q-operator in many cases (see, e.g. [35] and references therein). The models we shall discuss in 
this article are the generic generalized r^^^^ -model (where no Q-operator is known), superintegrable 
CPM, and XXZ spin chain of higher spin at the roots of unity. In the present paper, we obtain 
two main results related to these models. First, by constructing the fusion matrix from the explicit 
fused L-operator, we derive the fusion relations with the truncation identity for these models with- 
out relying on the theory of Baxter's Q-operator. Second, by the algebraic Bethe ansatz of XXZ 
spin chain of higher spin at the roots of unity and certain special r^^^-models which include the 
superintegrable CPM, we show the s/2-loop-algebra symmetry of the root-of-unity XXZ chain for 
a higher spin as in the spin-^ case [15, 18, 19]. Furthermore, they all share the same structure as 
the superintegrable CPM about the Bctlic equation and the evaluation (Drinfeld) polynomial for 
the symmetry algebra as in [41, 42]. As a consequence, a conjecture raised in [36] about the simple 
polynomial property for the s/2-loop-algebra evaluation parameters in the XXZ chain of spin-^^^^ 
has been justified. Note that for the r^^^ -model in the generic case where algebraic Bethe ansatz 
cannot be applied due to the lack of peudovacuum state, our fusion-matrix study strongly suggests, 
with computational evidences in cases, that the fusion relations with the truncation identity always 
hold. Hence by [24], the separation-of- variables method provides the solution of Baxter equation 
associated to the r^^^-model. 

The quantum inverse scattering method/algebraic Bethe ansatz developed by the Leningrad 
school in the early eighties [23, 28, 33] systematized earlier results about the Bethe ansatz of 
two-dimensional lattice models in an algebraic scheme by using the Yang-Baxter (YB) equation 
as a central role of solvability. A YB solution defines a local L-operator, which gives rise to the 
algebra of quantum monodromy matrices, called the ABCD-algebra. From the ABCD-algebra, one 
can construct a set of commuting transfer matrices, which in principal could be simultaneously 
diagonalized using a basis derived from the pseudovacuum state by the Bethe-ansatz technique. 
Furthermore, one can define the quantum determinant of the algebra, a concept first introduced in 
[26] (or see [28] Chapter VIII) , and played an important role in deriving the fused transfer matrices 
in this work. It is known that this algebraic method has long been used in the investigation of 
XXZ spin chains (see e.g., [31, 44] references therein), and in the r^^^-model [45, 46]. For the XXZ 
spin chain of higher spin in the root-of-unity case, it possesses some extra symmetry carrying the 
"evaluation" parameters, which indeed determine the eigenvalues of fusion matrices and implicitly 
encode symmetry of the model. Hence a new structure, not seen in the general XXZ chain of 
higher spin, appears in the root-of-unity theory. In this article, we employ the ABCD-algebra 
method in the generalized r^^^ -model and the root-of-unity XXZ chain of higher spin to study the 
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fusion matrix through some explicit fused L-operators. The boundary fusion relation will be our 
main concern. The technique is first to make use of the quantum determinant of L-operators, not 
only on the explicit form, but also its nature in commuting the fusion-product of elements so that 
the recursive fusion relation holds. Next the detailed analysis about "averaging" the L-operator 
leads to the boundary fusion relation. Since the ABCD-algebra of generalized r'-^^-model carries 
a non-equivalent, though similar, structure as the algebra for the XXZ spin chain due to the non- 
symmetric Boltzmann weights in the i?-matrix of r^^^-model, we shall provide a more elaborate 
discussion about the algebraic Bethe ansatz of the r'^^-' -model (though many like routine exercises in 
the field) . This is because the correct formulation with the precise expression of physical quantities 
is non-trivial, and required for the later CPM algebraic-Bethe-ansatz discussion in this work when 
comparing it with the complete results of superintegrable CPM derived from the functional relations 
[1, 4, 5, 10]. It is also needed for the parallel symmetry discussion between superintegrable iV-state 
CPM and the XXZ chain of spin-^^^^. The algebraic Bethe ansatz is known to be applied to the 
superintegrable r*^^^-model; however to what extent the results obtained by the algebraic-Bethe- 
ansatz method compared with the complete genvalues and its degeneracy known in the study 
of CPM [1, 5, 40] has not been fully discussed in the literature to the best of the author's knowledge, 
especially about possible symmetry structures of the model. To this end, we propose a scheme for 
certain special classes of generalized r^^^-model, with the superintegrable CPM included, where 
the pseudovacuum state exists so that the algebraic-Bethe-ansatz technique can be performed in 
the way like the root-of-unity XXZ spin chain. When applying to the superintegrable r^^^ -model, 
the setting enables us to conduct exact investigations of various problems. One can rediscover 
the Bethe equation, fusion relations, forms of eigenvalue spectrum and evaluation polynomials, 
known in the theory of CPM. Furthermore, certain eigenvectors derived from the pseudovacuum 
state can also be extracted by the algebraic Bethe ansatz method. Nevertheless, only certain 
sectors of the spectrum are covered by this scheme. In the case of the root-of-unity XXZ spin 
chain of higher spin, the algebraic Bethe ansatz method produces the correct form of evaluation 
polynomial for the degeneracy by a detailed analysis of the eigenvalues of fusion matrices. Then 
the zero-averages of off-diagonal elements in the quantum monodromy matrix, corresponding to the 
vanishing property of the AT-string creation operator, give rise to the sZ2-loop-algebra symmetry of 
the root-of-unity six vertex model by a "g-scaling" procedure in [15]. Thereupon one can identify 
the evaluation polynomial of sZ2-loop-algebra representation for a Bethe state through the explicit 
Fabricius-McCoy current ([21] (1.37)) of the model. 

This paper is organized as follows. In section 2, we discuss the fusion relations of the general- 
ized r^^^-model. We begin with some preparatory work in subsection 2.1 on the algebraic structure 
derived from YB relation for the generalized r^^^-model [45, 46]. Using standard techniques in 
the ABCD algebra and quantum determinant for the twisted i?-matrix, we construct in subsection 
2.2 the fusion operators from the fused L-operators so that the recursive fusion relation holds. 
By studying the average of L-operators, we then show evidences, verified in cases by direct com- 
putations, that the boundary fusion relation is valid for the generalized r^^^-model in subsection 
2.3. In section 3, we study two special classes of BBS models, which include the superintegrable 
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CPM, by the algebraic-Bethe-ansatz method where the pseudo-vacuum exists. We then perform 
the investigation on the Bethe equation and Bethe states for such models in subsection 3.1. The 
algcbraic-Bethe-ansatz discussion of special BBS models when restricted on the superintegrable 
7-(2)_]-Qodcl recovers the Bethe equation and evaluation polynomial of Onsager-algebra symmetry in 
the superintegrable CPM [1, 5, 40]. The comparison of those algebraic-Bethe-ansatz results with 
the complete results known in the theory of superintegrable CPM is given in subsection 3.2. In 
section 4, wc study the root-of-unity symmetry of XXZ spin chain with a higher spin. First we 
briefly review some basic concepts in the algebraic Bethe ansatz of XXZ spin chain that are needed 
for later discussions, (for more detailed information, see e.g., [23, 33] references therein). Then we 
summarize results in [27, 30, 32, 35, 42, 43] about the fusion relation for the spin-^ XXZ chain at 
roots of unity. Using the fused L-operators, we extend the construction of fusion operators in the 
spin-i to the spin-^=J^ XXZ chain at ATth root of unity for 2 < d < iV in subsection 4.1, where 
the fusion relations are derived. Furthermore through the fusion-matrix-eigenvalue discussion, we 
extract the correct form of evaluation polynomial incorporated with the Bethe equation. In sub- 
section 4.2, we show that the root-of-unity XXZ chain of spin-^^ possesses the .s/2-loop-algcbra 
symmetry, and verify the evaluation polynomial by the explicit Fabricius-McCoy current of Bethe 
states. In subsection 4.3, wc make the comparison between the root-of-unity XXZ chain of spin- 
and the A''-state superintegrable CPM, which are known to be closely related in literature 
[1, 5, 36]. We close in section 5 with some concluding remarks. 

2 Fusion Relations and Algebraic Bethe Ansatz of Generalized 
7-(2)_i2iodel 

We first briefly review some basic structures in the ABCD-algebra for the generalized r^^^-model 
in subsection 2.1. Then in subsection 2.2 we construct the fusion operators as the trace of fused 
L-operators so that the recursive fusion relation holds, and the boundary fusion relation will be 
discussed in subsection 2.3. 

2.1 ABCD-algebra and quantum determinant in the generalized r^^^^-model 

We start with some basic notions about algebraic structures in the generalized r^^^ -model. The 
summary will be sketchy, but also serve to establish notations, (for more detailed information, see 

[45] and references therein). 

For a positive integer A^, we fix the Nth root of unity, to = e n , Denote by C the vector 
space of A-cyclic vectors with {|^)}„gZjv standard basis where Z^v = Z/NZ , and X, Z the 

C'^-operators defined by X\n) = |n-|- 1), Z\n) = a;"|n) for n G Zjv, which satisfy the Weyl relation, 
XZ = uj~^ZX, with X^ = Z^ = 1. The L-operator of the generalized r^^^-model is built upon 
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the Weyl operators X, Z with C-auxihary space and C^-quantum space^ : 



L(i) 



1 + tKX (7 - 5X)Z 
t{a-(3X)Z-^ ta-i+^X 



m Bit) 

C{t) Dit) 



te C, 



where a, I3,j,S,k e C are parameters, which satisfy the YB equation^ 

i?(Vt')(L(i)(g)l)(l(g)L(tO) = (l(g)L(t'))(L(t)(g)l)i?(tA'), 



(2.1) 



(2.2) 



for the i?-matrix of the asymmetric six-vertex model, 

/ - 1 

t-l uj-l 
t{uj - 1) (t - 1)0; 



R{t) 



\ 





Then the monodromy matrix for the quantum chain of size L, 



(g)L^(i) =Li(i) Li (t) 

^=1 



AL(t) BL{t) 

cut) DL{t) 







tu-l 



L^{t) := h{t) at site 



(2.3) 



(2.4) 



again satisfy the YB equation (2.2), and the w-twisted trace, 

r(2)(t) := AL{0Jt) + DL{ut), 



form a family of commuting operators of the L-tcnsor space C of C . We shall denote the 
spin-shift operator of ® again by X{:= Y\i=i -^i) if confusion could arise, which carries the 
Z^r-charge, denoted by Q = 0, . . . , — 1. By 

[X,Al] = [X,Dl]=0, XBl = u;-^BlX, XCl = ^ClX, (2.5) 

X commutes with the r'^^-'-matrix. The relation (2.2) for the monodromy matrix gives rise to 
an algebra structure of the operator-entries A^it), BL{t),CLit), D^it), called the ABCD-algebra, 
(algebra of quantum monodromy matrix), in which the following conditions hold: 

[Ait), Ait')] = [Bit), Bit')] = [Cit),Cit')] = [Dit), Dit')] = 0; 

itiu - t')Ait)Bit') = it- t')Bit')Ait) + tiio - l)Ait')Bit), A,B — >C,D 

itu - t')Bit)Ait') = it- t')u}Ait')Bit) + t'((j - l)Bit')Ait), A,B — >C,D 

itu - t')Cit')Ait) = it- t')Ait)Cit') + t'iiv - l)Cit)Ait'), A,C — >B,D 

itu - t')Ait')Cit) = it- t')ujCit)Ait') + tiuj - l)Ait)Cit'), A,C — > B, D. 



(2.6) 



''Here we use the form of L-operator in accord with the convention used in [10, 40], which is essentially the 
transpose of the L-operator in [24, 45]. 

®Notc that (2.1) satisfy the YB relation (2.2) as well for a general u> not necessary a root of unity, using the Weyl 
operators X, Z with XZ = ui'^ZX. 
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The quantum determinant follows from the ABCD algebra by setting t' = ut: 

D{t); 

(2.7) 



B{ut)A{t) = A{ut)B{t), D{ujt)C{t) = C{ujt)D{t); 

A{t)C{ujt) = ujC{t)A{ujt), B{t)D{ut) = ujD{t)B{ujt), 

detq{® U){t) := D{ujt)A{t) - C{ut)B{t) = A{ujt)D{t) - B{ujt)C{t) 

= A{t)D{ut) - uC{t)B{ujt) = D{t)A{u;t) - oj-^ B {t)C {ut) . 



or equivalently, the quantum determinant of the monodromy matrix (2.4) is characterized by rank- 
one property of R{u)~^) in the following relation, 

R{u;-^){(S)Li{t)(^l){l(^(^Le{Lut)) = (1 (g) «)L£(a;i))(«)L^(i) (g) l)i?(a;-^) = detg(®L^)(i)-i2(a;-i), 

aux aux aux aux 

(2.8) 

with the explicit form for detq{(SiL^){t): 

detq{<^Li){t) = q{t)^X^, q{t) := — + {ad + ul3-f)t + uja-ynt^ . (2.9) 
The third- and fifth relations of (2.6) yield 

mB{s) = ^^Bis)A{t) + ^^Bit)A{s), Dit)B{s) = ^B{s)Dit) - ^^B{t)Dis). 
By moving -B(tj)'s to the left hand side of A{t), D{t), one obtains 



A{t) n-=i Biu) = n-=i 1^ ■ nr=i B{tMit) 



+ ET=i Mr nrii.^fc ■ Bit) nr=i.^. B{tMitk) 



D{t) a=i Biu) = nrii zth ■ n^i Bitmt) 



(2.10) 



T.T=i %^ nrii.^. 1^ ■ Bit) 1.^. B{tmtk) 



Similarly the second- and fourth relations in (2.6) yield 

Ait)ur=iCiti) =mlT^■uT=lCit^)Ait) 



T-rm t—ujtA T-rm. ^/i \r^/i\ \ ' J 



Dit)uZiCiti) =iiT=i'-f^-UT=iCiti)Dit) 

{U!-I)tk T-rm tk-Ultj 



+ Er=i ^ ni^i.^. ■ cit) nr=M^fc cit,)Dit 



■J) 



Note that by scaling the t-variable, parameters in (2.1) can be reduced to the case a+'j = 0, among 
which, with one more constraint u)(3 + 5 = 0, one can express —a = 7 = y~^, —(5 = oj~^5 = 
— K = fi^y^"^ for ix,y,fi) G C^. For the iV-state CPM, the rapidity variables of L-operator (2.1) 
are defined by 

kx^ = 1- k'fi~^, ky^ = 1- k'n^, {x, y, n) G C^ (2.12) 

where k' ,k arc temperature- like parameters with k'^ + k'"^ = 1. In the superintegrable case, the 
parameters in (2.1) and the quantum determinant are given by 

-a = -/3 = 7 = = -K = 1, detgLit) = Loh'^ it) X, (2.13) 

(see, e.g. [40] Sect. 5). Hereafter we shall always use h{t) to denote 

h{t):=l-t. (2.14) 
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2.2 Fused L-operator in generalized r^^^-model 

Here we construct the fused L-operator L^j\t) for the fusion r^-'^-matrix with L^'^\t) = L(t) in (2.1). 

For convenience of notations, we shall also denote the standard basis | ± 1) of the C^-auxiliary 
space of L{t), and its dual basis hy x = |1) , y = | — 1); x = (1| , y = (— 1|. For non-negative 
integers m, n, we denote by x"^y^ the completely symmetric (m + n)-tensor of defined by 

X y =x(8)...(8)x(8)y(8)...(8)y + all other terms by permutations, 
n I ' . ' 

^ m n 

similarly for x"*y". For j > 1, the C-'-auxiliary space is the space of completely symmetric (j — 1)- 
tensors of the C^, with the following canonical basis e^^^ and the dual basis e^-'^*: 



(i) _ ^j-i-k^k Jj)* _ ( J - 1 - 



X-' y , ej^ 



k 



X 



:^-^-V, fe = 0,...,i-l. (2.15) 



The L^-') [t) is the operator i (t) ) with the C-'-auxiliary and C^^-quantum space, where 

V ' J 0<k,l<j-l 



L^-'l(t) is expressed by 



LfejlO = {e'i^*U^^~^t) ®aux L(wt) ®aux ^{t)\ef^)- (2.16) 

Then lP\t) are intertwined by some i^^-'^-matrix. With lS-^\t) as the local operator, its monodromy 
matrix defines the commuting family of r^-'^-operators of ® C^, 

r(^)(t)=tr(..((g)Lf (a;i)). (2.17) 

1=1 

We now show the fusion relation between t^~^^\ t^^^ and t^^~^^ through the quantum determinant 
(2.9). 

i+i 

Consider the auxiliary-space tensor (g) as a subspace of (g) with the identification 

and denote .= x (g) e)^^-^ - y ® ej: for < < j — 2. Then % form a basis of 

(g) C-^ with the dual basis ep"'"^''*, /j^"'"^''* expressed by 

Then the expression of , el yields 

^4.+i)*|LO+i)(,)|eO+^)) = {e^^+'>\L{u^-H) 0_ LO)(i)|ep+^V (2-18) 
In order to determine the rest entries of L{LU^~^t) (Siqux ^^■'\t), we need the following simple lemma. 
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Lemma 2.1 The second equality in (2.8) is equivalent to the following relations: 

(x2|L(a;t)Oa«a;L(t)|xAy) = {y'^\L{ujt) ®auxl^{t)\x hy) =0, 

(x (g) y|L(wi) (g)a„a; L(t)|x A = {y ® x\h{tjjt) ®aux'^{t)\ - x ^y) =\deiqL{t), 



(2.19) 



(2.20) 



where x Ay = ^(x^y — y(8)2). Hence we have 

{e^^'^*\L{uot)®auxUt)\x®y) = {e'j^^*\LiLjt)(E)auxHt)\y(S)x) for fc = 0,l,2. 
As a consequence, for an integer j > 2, and Vi = xory for 1 < i < j — 1, we have 

{e''^^*\L{iO^''^t) ®aux ■ ■ ■ ®aux L(wt) ®aux L(t)|fl ® V2 ® ■ ■ ■ ® Vj-i) 
= {el^^*\L{uJ^~'^t) <Siaux ■ ■ ■ ®aux L(wi) ®aux L(i)l^<Ti ® Vcr^ ® ' ' ' ® Vaj_-^) 

for < k < j — 1, and all permutations a. 
□ 

Since the entries of L{u}^^^t) 0aux l-^^\t) are determined by those of L{ui^~^t) 0aux ■ ■ ■ ^aux L(i), by 
(2.20) one has 

(egi')|L(a;^-4) ^aux l^^'\t)\fj^?) =0, -1 < k < j - 1, < k' < j - 2. (2.21) 

Using (2.19) and (2.20), one finds 

{fj^'-'^*\L{co^-H) ®aux L(^-)(t)|/ir')) = {fl^-^^*\L{oj^-H) ®aux L^^\t)\{2x A y) y>^') 

= T^iii'l')^^ <^y- ii+\)y ® 2:) ei^''>\L{u;^-H) ®aux l^^'\t)\{2x Ay)® x^-'^'-^ ® y^') 

\k+l) 

= det,L(a;^-2t) • {e'j^-^>\Li^-^){t)\x^-'''-^ ® y^'), 
which, by (2.9) and (2.19), in turn yields 

{ft'>\Liu;^-H) ®_ L(^)(t)|/(r^)) = {et'^*\L(^-'Ht)\e^^r%{u^-'t)^ (2.22) 

Prom the definition (2.17) of r^-'^-matrices, the relations, (2.18) (2.22) and (2.21), imply the follow- 
ing result: 

Proposition 2.1 The t^^^ -matrices satisfy the (recursive) fusion relation by setting t^^^ = 0, t^^^ = 
I, 

r(2)(a;J-it)rO)(t) = ziuj^-H)XT(^-^\t) + T^^+^^t), j > 1. (2.23) 
where z{t) = q{t)^ with q{t) in (2.9). 
□ 
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2.3 Boundary fusion relation in generalized r^^^-model 



For convenience, we introduce the following convention for a family of commuting operators 0{t): 



n-l 



[0]n{t) := n 0{ujH), n e Z>o, 



i=0 



and the average of 0{t) is defined by 

(o) (= (o)(t^)) = puit). 

The "classical" L-operator of BBS model is the average of (2.1): 



.N 



^ ^ '(C) {D) I (-l)^+i(«'^-/3^)^^ ^ + (_l)iv+i^iv^iv^iv I- ^ - ) 



The averages (^l), (-Bl), {Cl), {Dl) of the monodromy matrix (2.4) coincide with the classical Lth 
monodromy associated to (2.24) [45]^: 

By (2.20), the [k, l)th. entry L^"'](t) of L(-'^(t) is equal to the expression in (2.20) with Vi = x for 
l<i*<i~l~^5 and y otherwise. Hence one can express L^^\t) in terms of entries in (2.1). For 
example, the matrix-form of l(^)(= L^^\t)) is 

/ [AUt), A{ut)B{t), [BUt) \ 

A{iut)C{t) + C{ujt)A{t), A{ujt)D{t) +C{ujt)B{t), D{ut)B{t) + B{ujt)D{t) . (2.26) 
V [cut), C{u;t)D{t), [DUt) ) 

Note that one can also write lJ^J = A{Lot)B{t),L'^} = D{ujt)A{t) + B{ut)C{t), L^^J = D{ut)C{t). 
Among the L^^^-entries for a general j, the following ones can be derived by setting ®Vi = x^~^~^®v^ 
in (2.20): 

l!il = [A],.,{t), Lg_i = [S],_i(t), l(^,,o = [C'b-iW, l^\,., = \D\,-i{t); 
L\^} = [A]^_i_^{Jt)[B]i{t), l!^^\^, = [C],.i.,{Jt)[D]i{t), l<l<j-2. 

By the relation, e): = e):_-y ® y + ® x, for the dual basis of C-^'s, one can compute L]^^ 
in terms of L^-'^-entries using the following recursive relations: 

= Lif2i,oM)^^W + Lg(cut)^(t), 



®The formula (2.25) about averages of the monodromy matrix is stated in [45] page 966 as a consequence of 

Proposition 5 (ii) there, (or Lemma 1.5 in [46]), but without proof, and also not with the required form which could 
be misprint in both papers. The correct version should be < A(r) >=< Ti >< T2 >, instead of A(< T >) =< 
Ti >< T2 >. As the author could not find a proof in literature about the correct Tarasov's statement, here in 
this paper we provide a mathematical justification about the correct statement in Proposition 2.2 and 4.1 for the 
generalized r'^^-model and the XXZ spin chain respectively. 



10 



by which, one can in principal derive the expression of L^-'](t) from those in (2.27) for a given j. 
However it is still not easy to obtain a close form of L^-';^ for all j except the following {k, Z)'s: 



(2.29) 



Indeed, the above formulas are derived by the induction-method using the relations, [C]k{i^t)A{t) = 
u;''A{iu-''t)[C]k{t), [B]k{ujt)D{t) = Lu''D{uj-H)[B]kit) for k>0, which follows from the equalities 

C{t)A{t) = A{uj-H)C{t), B{t)D{t) = ujD{u-H)B{t). (2.30) 

In the study of CPM [10], the r^^^-fusion relations (2.23) with rapidity parameters in (2.12) are 
truncated at t^^'^^'> with the following boundary fusion relation: 

r(~+i)(t) = z{t)XT^^-^){ut) + u{t)I, u{t) := {{Al) + {DlW), (2.31) 

where Al,Dl are operators in (2.4). We are going to indicate the relation (2.31) always holds 
in the generalized r^^^-model for arbitrary parameters. Express the diagonal entries of (L + l)th 
monodromy matrix in (2.4) by ^L+i(t) = (tIl A + C)(t), BL+i{t) = {Al B + Bl D){t) 
, {Al+i,Bl+i,Al,Bl Cl+i,Dl+i,Cl,Dl); hence their averages: {Al+i) = [Al <^ A + Bl <^ 
C]N{t), {Bl+i) = [Al ^ B + Bl D)]N{t), etc. Since Al,Bl and Cl,Dl satisfy (2.7), using the 
expression (2.16) for L^.^'*'^'' , L^,^^-* , one obtains 

{Al+i) = {Al) {A) + {Bl) (C) + Y.k=ML]N-kiioh)[BLUt) 4^0^'^ (t), 

(Bl+i) = (Al) {B) + (Bl) ® {D) + Ef=Y[^L]iV-fe(^'=*)[5L]feW ® 4!!?'^ (2-32) 

Al+1, Bl+1, Al, Bl — > Cl+i,Dl+i,Cl,Dl. 

Using the above relations and the following proposition on the vanishing L^^"*"^^ , L^.^^-* for k ^0,N, 
one can determine the average {Al), (Bl) etc. through the relation (2.25). 

Proposition 2.2 The entries Lj,^^^-*, L^^"*^^ ofL^^'^^\t) (on the 1st and Nth column) are all zeros 
except 4^"^^^ = {A), L^^^^^ = [C), Lq^^^ = [B), L^^^^^ = [D). As a consequence, the relation 
(2.25) holds; hence u{t) in (2.31) is equal to tr(>C(t^)^) with C(t^) in (2.24). 

Proof. We need only to show the vanishing of L^^"*"^-*, L^^^-* for k ^0,N. The relation (2.32) with 
L = 1 yields 

(A,) = {Af + {B){C) + E^=M]N-k{u;h)[BUt) ^^^^(t), 
(B^) = {A){B) + {B){D) + E^-MN-k{<^''t)[BUt)®Li';j^'\t), 

which are operators invariant when changing t by cot. Note that X'^Z^ (i, € Zjv) form a C-basis 
of all A^-by-A^ matrices; and the term [A\N^k{'^^t)[B]k{ujH) in the above expressions can be written 
in the form J2m=o Pmit)X™' for some Pm{t) £ By this, each term appeared in the above 

expression of (^2), (-S2) is again invariant under t 1— > ut, i.e., for 1 < A; < A?^ — 1 and e = 0,N, the 
equality holds: 

[AU-k{u:''+H)[BUu;t) ® 4^/'^(a;i) = [AU-kiu;''t)[BUt) ® 4f,+'^(0- 
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This implies ^(t)(8)4^+^^(a;t) = A{u;H)®L[^^'^^\t) by the relation, A{coH)[A]N-k{i^''^^t)[B]k{iot) = 
A{t)[A]N-k{u>^t)[B]k{t). The linear independence of matrices 1 and X in the expression of A{t) in 
turn yields 

Therefore L[^o^^\t) = L[^N^\t) = fov 1 < k < N - 1. □ 

Remark. For the superintegrable r'^^^-model (2.13), by (2.24), the above proposition implies the 

averages of Lth monodromy matrix-entries are: {Al) = (Dl) = (1 — t^)^, (Bl) = {Cl) = 0. In 
CPM where the rapidities p,p' of r^^^ -model are in (2.12) with p fixed and t = Xp/ypi for p', one 
can compute the eigenvalues of (2.24) using the coordinate A = /j,^ to express the u{t) in (2.31) for 
CPM: uit) = ap{Xp,)+ap{XpO where apiXp>) = (^^i^^^)^ ( [10] (4.27)). 

By Proposition 2.2, the {N + l)th fused L-operator can be written in the form 



■ (N+l) 



it) 



( * {B) \ 

L(^-i)(t) 



where * = (4^^^ • • • ,4^/iJ), *' = (4^,+^ • • ■ ,L^^^^]},), K^-^Kt) = (4^^^ W) with 

L^^l~^\t) := L^^^J^^(t). The boundary fusion relation (2.31) is equivalent to the relation between 
7" (^-1) -matrix and the Lth monodromy matrix for L^^^^^: 



L 

z{t)XT^^-^\Lot)=trfjN-,i(^Lf-^\u;t)). (2.33) 



=1 



The above relation (2.33) will hold if we can prove L^^-^\ujt) similar to L^^-^\ujH) ■ detgL(t) by 
a diagonal similarity transformation M (independent to t): 

(ut) = Ml(^-^) {ujH)M-^ ■ detgL(t). (2.34) 

We now determine the above similarity relation for N = 2,3,4. For N = 2 where u = —1, by 
(2.26) one has L^^^ut) = A{t)D{ut) - ujC{t)B{ujt) = detgL(t), hence follows (2.33). For AT = 3, we 
use (2.20) to express L'^-^+^^-entries, 

L{,2)(a;t) = A{t)A{uj^t)D{Lot) + A(t)C{ujH)B{ujt) + C{t)A{LoH)B{ujt), 
L(^j(u;t) = B{t)B{ujH)C{ujt) + B{t)D{ujH)A{ujt) + D{t)B{ujH)A{ujt), 
ifliut) = A{t)C{ujH)D{ujt) + C{t)A{ujH)D{ujt) + C{t)C{ujH)B{ujt), 
l}il{ujt) = B(t)D{ujH)C{ujt) + D{t)B{ioH)C{ujt) + D{t)D{ioH)A{Lot). 

By (2.30) and (2.9), one arrives 

L[)^)(u;t) = A{ujH)A{t)D{ujt) + {LoA{t) + uj'^A{ujt))LoC{t)B{Lot) = A{ujH)detgL{t), 

Eg (cut) = B{uH)B{t)C{ujt) + B{t){D{ioH) + lj^ D{Lot))A{Lot) = -u;B{ioH)detqL{t), 

Lf)^{ujt) = C{ujH){A{u;t) + ujA{ujH))D{ujt) + C{u;H)C{t)B{ujt) = -uj-^C{LoH)detqL{t), 

L^liujt) = {uD{ut) + D{t))B{t)C{ut) + D{ojH)D{t)A{ut) = D{ioH)detgL{t). 
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Hence (2.34) holds with M = dia[— a;, 1]. Indeed, the same method as in above arguments, together 
with (2.27) and (2.29), yields the following four relations for an arbitrary N with a;-^ = 1 : 



For = 4 where lu = using the explicit form of A,B,C,D in (2.1) to compute the 

L{3)_cntrics other than those in (2.35), then comparing L^^^ with (2.26), we can verify the relation 
(2.34) holds with M = dia[— \/— Ij ^~~2^' ■'■^ arbitrary given N, the similar transformation 

in (2.34) could be obtained by direct calculation, however the general structure has not been found 
yet. 

3 Algebraic Bethe Ansatz in Generalized r^^^-model 

In this section, we use the algebraic-Bethe-ansatz techniques to discuss certain BBS models centered 
at the superintegrable r^^^^-model, where the pseudovacuum state exists. 

First we note that the non-trivial kernel of the operator C{t) in (2.1) is equivalent to (C) = 0. 
By scaling the t-variable and changing X by uj'^X, one may set a = /3 = —1, i.e.. 



Alt) / 1 + i-,-,X)Z 

' C{t) D{t) J \-t{l-X)Z-^ -tj-^X ' ^ ' 



In this section we shall consider only the case (3.1) with the following constraint: 

7 / u;-''6 for /c = 2, . . . , A^, (3.2) 

which will be called the special BBS model in this paper. The quantum determinant in (2.9) now 
takes the form 

detqL{t) = q{t)X, q{t) = Lohi{t)h2{t), where hi{t) := {I + Kt), h2{t) := {-jt ). (3.3) 

AC 

The special BBS model with one further constraint (A) = (D), equivalently 7^ = = {—k)^, 
implies 7 = u;~^5 = —lu~'^°k for some zq. By substituting 'yt by t, then applying a similar trans- 
formation, L(t) is equivalent to the case 7 = 1, with B{t) = (1 — ujX)Z, A(t) = 1 — tu'^^'tX, 
D{t) = -t + ui^°X, and hi{t) = ui^'^h2{t) = h{t), where h{t) is in (2.14). In particular, when io = 0, 
one arrives the superintegrable r^^^ -model (2.13). 

3.1 Bethe equation of the special BBS Model 

Denote the eigenvectors of the C^-operator X by 

fk := Z\ Yl I"))' ^(/fe) = '^"'/fe for k G Zat. 

Let Al, Bl, Cl, Dl be the entries of the Lth monodromy matrix (2.4) with L(t) in (3.1) satisfying 
the condition (3.2), and the transfer matrix is given by T^'^\io~^t) = AL(t) + DL{t). It is easy to see 
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the C{t) in (3.1) has one-dimensional kernel space generated by v := fi, and A{t)v = h\{oj ^t)v, 
D{t) = h2{t)v, with hi,h2 in (3.3). Define the pseudo-vacuum: 

O.L = 0i=ive C^, ve = v for all i, 

then AL{t)nL = /n(a;-H)^OL , DL{t)9.L = h2{t)^^L, CL{t)VlL = and 

BL{t)^L = (7 - <^w^) Eii hi{u-HY-^h2{t)^-^{l ®---®Ze,th®®---® l)^^- (3.4) 

Note that for L = 1, we have = fi and BL{t) = B = - 5X)Z. Then B^^l = Y{^i=2{l - 
Slo'''^~^) fk+i for < k < N — 1 are non-zero vectors by (3.2), which form a basis of C^. 
For m distinct non-zero complex numbers tj,l < j < m, we consider the vector 

fn 
i=l 

and define the t-function, 

A(t;ti, . . . ,t^) = h.iu-H)'^^^^^ + h2{t)'^u-^^, Fit) := ^(1 - tTH). (3.5) 
The polynomial-criterion of A{t;ti, . . . ,tm) is given by the following Bethe equation for tj's: 

— — w — r^, z = l, ...,m. (3.6) 



By (2.5), the Zjv-charge of pseudo- vacuum and Bethe vectors are 

XnL=u;-''nL, X^{ti,...,tm)=u;-''-'^^{ti,...,tm). (3.7) 

Proposition 3.1 The polynomial condition of A{t;ti, . . . ,tm) is equivalent to the Bethe equation 
(3.6) for ti 's, which is the necessary and sufficient condition for ^(ti, . . . , tm) as a common eigen- 
vector of the transfer matrices r^^^ (^^"^0 ■ In this situation, A(t; ti, . . . , tm) is the eig en-polynomial 
for the eigenvector ^'(ti, . . . , tm)- 

Proof. Write the function (3.5) in the form 

Kit- ti, . . . , tm) = h,{u-H)^ n ferr + ^2^*)" fi i^^^y 

By (2.10), {AL{t) + DL{t))]\T=i BL{ti)QL is equal to the sum of A(i;ti, . . . ,t„,)*(ii, . . . ,tm) and 
Afe(i; ti, . . . , tm)B{t) Y{T=i,i^k BL{ti)nL, ik = l,...,m) where 

^*'--->*")-S(M-'*.)^n|f^-M«'n^). 

By this, ^'(ti, . . . , tm) is a common eigenvector of t) provided that Afc(t; ti, . . . ,tm ) = for 

all k, which is equivalent to the Bethe equation (3.6). Then follows the proposition. □ 
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Remark In the above algebraic-Bethe-ansatz discussion of BBS model, we set C{t) = —t{l — 

X)Z^^, which is the one appeared in superintegrable CPM. One may also set B{t) = (1 — ojX)Z 
instead, and let C{t) = t{a — f3X)Z~'^ with parameters a, /3; then conduct a similar algebraic- 
Bethe-ansatz discussion. In the latter setting, the pseudo- vacuum Q'^ = 'S>{fo)e satisfies the 
relations, BLit)^'^ = 0, ^lW^^l = hi{t)^n'^, and DL{t)n'^ = ojh2{u-H)^l'^. With the same 
argument in Proposition 3.1 but using the relation (2.11) instead, one considers the Bethe vector 
...,tm)-- rii^i CL(ti)n'^, with the Bethe equation 

which in turn yields the eigenvalue expression for ^'^{ti, . . . ,tm)- 

A\t;t„ ...,tm) = u:^^-h,{u-Hf^^^^ + hi{t)'^. (3.9) 

□ 

We now consider the special BBS model (3.1) in AT = 2 case where u = —1, and the Bethe 
equation (3.6) becomes 

(^:^^)^ = (-l)™+\ z = l,...,m. 

The solutions of the above relation are determined by non-zero roots of (7^ — 6/ k)^ = ±{Kt — 1)^, 
each has L solutions when ^ ±k^ ^ 6^: 

fk = —, ru = —, — , ^^i,r % s [1 < k < L), gl e l . 

k(7— it2k) k(7— (— 1)-^/-^(t2k) V — — /' 

L „ 

For generic 7, k and S, the Bethe states are expected to generate the quantum space (g) C , e.g., 
for L = 2, it follows from the expression of the Bethe states Cl2{= /i <8) /i), ^'(ri), ^'(r2), '^{r[,r2): 

^^iri)= hi{-n)h®fo + h2{ri)h®h, (i = l,2); 
(:F5F*W,^2) = ^(/ii(-r-i)/ii(-r^) + /i2(r'i)/i2(r^))02 + (/ii(ri)/i2(r^) - /i2(-rl)/ii(-r^))/o ® fo 

Note that in the superintegrable case (2.13), is the only Bcthc state in the above setting. Indeed 
in this situation, T^'^\-t) = (1 + i^)(l + X (g) X) - 2t{Z (g) Z){1 ~ X (g) X), with the eigenvalues 
2(1 + t^), ±4t, and cigcnspaccs (^^2, fo fo), {fo "X" /i T /i "X" fo) respectively 

For > 3, it is not easy to obtain the complete solutions of the Bcthc equation (3.6) in general, 
hence difficult to determine the space generated by Bethe states in (g) C . For m = 1, the Bcthc 
equation (3.6) becomes h2{t)^ = hi{LJ~^t)^ , which has L solutions in generic cases. By (3.4), the 
Bethe states for m = 1 likely form a basis for the subspace generated by (1(g)- • -(gZith®^ ■ ■ ■^1)0,'" 
for all i. 
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3.2 Algebraic Bethe ansatz in superintegrable chiral Potts model 

We now discuss the superintegrable r*^^) -model (2.13) in the setting of previous subsection, then 
compare the result with those in the theory of superintegrable CPM [1, 5]. 

In the study of the chiral-Potts-transfer matrix in superintegrable CPM, there are quantum 
numbers. Pa and Pj,, appeared in the eigenvalues, which satisfy the conditions 

0<r {■.= Pa + Pb)<N-l, Pb-Pa = Q + L (mod TV), LPb = m{Q - 2Pb - m) (mod iV) 

with Q the ZTv-charge as before ((C.3) (C.4) in [1], (6.16) in [5], or (59) (63) in [40]). The Bethe 
equation of superintegrable CPM is given by ((4.4) in [1], (6.22) in [5]): 

Hcp{uJ-hi)^ _ _, F{aj-Hi 



^ ' ^Z {i = \,...,m), F{t):=W{\-t-H), (3.10) 



where Hcp{t) is the polynomial defined by 

l-t^ 1 _ t^v 
«CP(*) = — (= (3.11) 

Define the normalized r^^)- and r^-^^-matrices ([42] (3.4)) by setting T(°) = 0, T^^\t) = Hcp{lo-H)^, 
and 

Then r(j)-eigenpolynomials determined by F{t) in (3.10) satisfy the relations ([42] (3.8)-(3.11)): 

T(^\t)F{t) = io-^Hcp{t)^Fico-H) + Hcpiio-H)^F{u;t), 

> IV 

(3.13) 



rO-)(t) = F{u-H)F{u^-H) ElJo %tiVA:^V ' ^ 



r(j)(t)r(2)(a;^-if) = u-^H{uj^-H)^T(^-^\t) + H{u^-H)^T(^+^^{t), {j > 1), 
r(^+i)(f) = a;-^r(^-i)(a;i) + 2H{uj-H)^. 

Note that the above T^-'^-relations are equivalent to the fusion relations (2.23) (2.31) for t^^^'s in 
the superintegrable r^^^-model ((53) in [40]). The second relation of (3.13) in turn yields 

N-l 

T^^){t) = u''''F{u-H)F{lo''-H) ^ Vcpit'') 

k=0 

where 'Pcp{t^) is the evaluation polynomial of the Onsager-algebra representation for the degen- 
erate eigenspaces in the superintegrable CPM 



Bethe equation (3.6) in the superintegrable r'^^^-model is the same as (3.10), but with the constraint 
r = Pa + Pfc = —m (mod N). The eigenvalue (3.5) of the Bethe state *(ii, . . . ,tm) associated to 
the Bethe solution tj's of (3.6) is characterized by the relation 



16 



The above relation is the same as (3.13) where the factor uj~^^ in formula (3.12) is identified with 
a;"*. By this, we arrive the conclusion that the r^^^-eigenspace for the Bethe state '^{ti, . . . ,tm) 
satisfies the constrain —m = Pa + = Pf, (mod N), hence lies in sectors with Pa = 0, —m = P^, 
hence —Q = L + m (mod A^). In particular, the pseudo- vacuum fi^ is the ground-state for 
= = with Q = -L (mod iV). 

We may discuss the algebraic Bethe ansatz of supcrintegrable CPM in the setting described in 
the remark of Proposition 3.1 with the pseudo-vacuum $7'^ and Bethe state '^'{ti, . . . ,tm) ■ The 
Bethe equation (3.8) now becomes (3.6) with r = —{L + m) (mod N). The eigenvalue (3.9) for 
'^'i^iti, . . . , tm) is expressed by 

Then the r^^^-eigenspace generated by ^'{ti, . . . ,tm) has the constraint, P5 = 0, —Pa = L + m 
(mod A^), hence Q = m (mod N). Now the pseudo- vacuum fi'^ is the ground-state for P^ = 
0, Pa = -L (mod N) with Q = 0. 

The above discussion about the superintegrable CPM has indicated that the genstates in 

the algebraic-Bethe-ansatz approach appear only in the sectors with P^ = or P^ = 0, not for the 
rest sectors. The full spectrum of r^^^ -eigenvalues and the symmetry of its degeneracy has been 
well studied in [40] through the chiral Potts transfer matrix as the Baxter's Q-operator for the 
r^^^-matrix, and the complete detailed structure was obtained by the functional-relation-method. 
Nevertheless, the algebraic-Bethe-ansatz approach does provide a mechanism to understand some 
-7-(2)_ eigenvectors in certain sectors. 

4 The XXZ Spin Chain of Higher Spin at Roots of Unity 

In this section we study the fusion relation and Bethe ansatz of XXZ chain of spin at roots of 
unity for positive integers d >2 with the anisotropy parameter q. Here the discussion of the XXZ 
spin chain will always assume with the even chain-size L. We shall denote u = q^; when q is a Nth 
root of unity, we consider only the odd N case, then take q, oj{:= q^), q^ all to be primitive Nth. 
roots of unity. 

The L-operator of the XXZ chain of spin-i is the operator with C^-auxiliary and C^-quantum 
space: 

Lis) = ( ) (.), . e C, (4.1) 

where the C^- (quantum-space) operator-entries Lj^j are given by 

-^^0,0 = 






with a, 6, c the g-dependent s-functions 

a = a{s) = sq~ — s~'^ q^ ^ b = b{s){= a{qs)) = sq'^ — s~^q~ , c=q — q~^. (4.2) 
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The operator (4.1) satisfies the YB relation, 

RUs/s')iLis)(^l){l(^L{s')) = {l(^L{s')){L{s)(^l)Reyis/s'). 

where i?6v is the i?-inatrix 



(4.3) 



aux aux 



( s-^q - .sg"^ 



i?6v(s) = 



v 









s q 
1 s- 




q — q ^ s ^ — s 







s^^q — sq~^ j 



Hence the monodromy matrix of size L, 

L 



e=i 



A(s) Bis) 
C{s) D{s) 



(4.4) 



again satisfies the YB relation (4.3), i.e. the entries A, B, C, D form the well known ABCD-algebra, 
in which the following relations hold: 

[A(.),A(.')] = ms)Ms')] = [C{s),Cis')] = [D(5),D(.')] =0; 

A{s)B{s') = fs,s'B{s')A{s)-gs,s'B{s)Ais'), A ^ B; A, B C, D; (4.5) 

D{s)B{s') = /,,,,B(s')D(s) - gs>,sB{s)D{s'), B ^ D; B, D 



A,C. 



where the functions fg^s' > gs,s' are defined by 



ss'jq^ - 1) 
q{s^ — s'^) 



(4.6) 



Then follow the relations 



A{s) nrii Bis,) = (n™ 1 fsMUT=i b(s.))a(5) - et=i 9s,k{nT=l,^^k a,)b(5) mu^^k Bis,)Ais,), 
Dis) uzi 6(^0 = (n™ 1 ks){UT=i b(s.))d(s) - Er=i 9kAmi,m Afe)B(s) 6(^00(5^); 

A(s) ^ D(s),B(s) ^ C(s). 

(4.7) 

Here we write only i for Sj in the subscripts of fs,s'igs,s'- 

As the matrix Rev{q) is of rank 1, the quantum determinant of the YB solution L{s) in (4.3) is 
defined by 

Rey{q){L{qs)<^l){l<^L{s)) = (1 L(s))(L(gs) (g) l)i?6v(g) =: detgL{s) ■ Rey{q), (4.8) 
or equivalently, 



aux aux 



Aiqs)Cis) = Ciqs)Ais), B{qs)D{s) = Diqs)B{s), 
B{s)A{qs) = A{s)B{qs), C(s)D(gs) = D{s)Ciqs); 
detgL(s) = Aiqs)D{s) - C{qs)B{s) = D(^s)A(s) - B(^s)C(s) 
= D{s)A{qs) - C{s)B{qs) = A{s)D{qs) - B{s)C{qs). 



(4.9) 
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For the local L-operator (4.1), detqL{s) = a{s)a{q^s), so the quantum determinant of (4.4) is 
equal to a{s)^ a{q^ s)^ . The statements in Lemma 2.1 are valid for the operator L{s) by replacing 
L(a;t) ®aux L(t),detgL(t) by L{s) ®aux L{qs),detqL{s) in (2.19), and L{iL!-^^^t) (^^aux ■ ■ ■ ®aux L(i) by 
L{s) ®aux ■ ■ ■ ®aux L{q^~'^s) in (2.20). Now we consider the transfer matrix of the six- vertex model 
of even size L, defined by the trace of (4.4). Denote 

r(2)(t) = s2i(A(s) + D(s)), t-qs"^ G C. 

Then T^'^\t) form a family of t-polynomial operators acting on the quantum space ® C^. As in 
the discussion of subsection 2.2, we constructed in [42] the fusion matrices T^^\t) (0 < j G Z) 
from the fused L-operator L^^\s) associated to L{s) in (4.1), which is a matrix with C^-quantum 
and C-'-auxiliary space, as follows. With the basis in (2.15) for the C'-auxiliary space, L^^\s) = 
(L^''](s) J where the C^-(quantum-space) operator L^^\{s) is defined by^ 

tU)^ \ _ (Sfc ^aux L{qs) (giaux ' ' ' ®aux L{q^ '"^ s)\ef^ ) 

U.i=o Kq's) 

The relations, (2.16) (2.21) (2.22), now turn to 

(ep+i)*|LO-)(,) L{q^-^s)\ft'^) = 0; (4.11) 
(/^^^1lO)(.) L{q^-^s)\ft'^) = b{q^-^s){et'^*\LU-^\s)\e^-'^). 

Prom the first relation in above, one obtains the explicit form of L^^\{s) for j > 2 by induction 
argument ([42] Proposition 4.2): L^'f!\{= L^i\{s)) {0 < k,l < j — 1) are zeros except k — I = 0, ±1, 
in which cases 

(,)_/a(g'=.) \ (,) _/0 q^-^-'' - q-^+^+>^ \ _/ 0\ 

'^'^ - 1^ a(g^-i-^s) ) ' ^-^''^ " V ; ' ^-^''^ " - r ' J 

The fusion matrix T^i\t) is defined by T(°) = 0, T^^\t) = (1 - u-H)^ , and for j > 2, 

r(^)(t) = {sq'^)Hvaux{<klLf{^))^ t = qs^e C. (4.12) 

The relations in (4.11) and the explicit form of L^"']'s in turn yield the fusion relation of T^^^'s with 
the truncation identity: 

r0)(i)r(2)(a;^-it) = (1 - u;^-H)^T(-^-^\t) + (1 - u^-H)^T(^+^\t), j > 1, 

T(^+i)(i) = r(^-i)(a;t) + 2(1 - u;-H)\ ^ ' ^ 

parallel to the fusion relations, (2.23) and (2.31), in the generalized r^^^-model. 



''One may use //'■''-operator defined by revising the order of L{s) ^aux L{qs) ®aux ■ • • ^aux L{(f~^s) in its definition, 
replaced by L{q^~^s) (^aux L{q'~^s) ®aux ■ ■ ■ ®aux L{s) as the order in (2.16). The latter form is related to the first 
expression in (4.8), instead of the second one in the former case. 
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4.1 Fusion relation and algebraic Bethe ansatz in XXZ spin chain of higher spin 



When interchanging the auxiliary and quantum spaces of the dth fused L-operator L^^\s) for a 
positive integer d, we arrive the L-operator of XXZ chain of spin-^=i, which is the matrix with the 
C^-auxihary and C'-quantum space: 

L(.)=f [''']is), sec, (4.14) 

V Li^o Ll,l / 

where the entries i\-a,h)i (0 < < d — 1) of La^b are zeros except (Lo,o)i = (Lia)^!^!} = a{q^s), 

(0 < i < d - 1), and (Lo,i)*+i = q'^-^-' - q'd+-^+\ (Li,o)*+^ = - q-'-\ (0 < i < (i - 1). Here 

\-a^h ai'e operators of the quantum space with the standard basis e* (i = 0, . . . , d — 1). It is 

well-known that the above expression of gives the c?-dimensional irreducible representation of 
Uqish) (see, e.g., [31, 32]), 

d-2 -1^ _, -(d-2) -^1 _ 

\-Q^Q = sq 2 K 2 - s 2 K2, Lo,i = (g - g ^)e_, 

Li,o = (g-g )e+, Li,i = 2 7^2 - s 2 2 , 



with the relations, Ke±K~^ = q^'^e± and [e+,e_] = ^gZf-i^ ; where if^ = g'S''' and the spin- 
operator of C^: = dia[^^, . . . , ~^"^^ ]. By the direct verification, one finds the L-operator 
(4.14) again satisfies the six- vertex YB relation (4.3), hence the Lth monodromy matrix 

(g) leis) = Li(s) ® • • • ® U(s) = ( ^[''^ ^^'l ] , le{t) := l{t) at site £, (4.16) 

f=i V c(s) D(s) y 

has the entries forming the ABCD-algebra with the pseudo-vacuum and quantum determinant 
(4.9): 

n =1 e° C'^', detgLL(s) = 0(9^^5)^0(5)^ • id. (4.17) 

The trace of (4.16) defines a family of commuting ((8) C )-operators, called the transfer matrix of 
the XXZ chain of spin-^^: 

t{s) = A{s) + D{s), 

which commute with q^'' =(8) where is the spin-operator of (8) C^. Similar to the discussion 
of (2.16) and (4.10), we consider the fused L-operator associated with (4.14)^: 

^k}(s) = (efc^*l'-(s) ^aux \-iqs) ®aux ' ' ' ®aux \-{q^-'^ s)\e^^) , 
and define the fusion matrix 

tO)(5) = tr(.,((g)Lf(.)), 



"^Note that for d > 2, l2\[s) have no common factors in contrast to the case d = 2 in (4.10). Hence we use 
the notation L for the L-operator for a general d to distinguish it from the L-operator (4.10) for d = 2, in which 
case \J'^\{s) = (^{Zt^^'i^ ^))^^k\^^)- However the fusion matrix t^^'' will still be divisible by certain factors after 
multiplying the normalization factors in (4.19) and (4.27). 
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with t^'^\s) = t(s). Then t^^\s) form a family of commuting operators of (g) C^. Similar to (4.11), 
we now have 

^^^^ L(g^-is)|//^-^)) = det,L(g^-2,)(ei;'-^)*|L0-i)(s)|ep-^)), 

then follows the fusion relation of the XXZ chain of spin-^^, parallel to (2.23) in the generalized 
T^^) -model, by setting t(°) = 0,t(^) = /: 

t(j)(s)t(2)(gi-is) = a(g'^+J-2s)^a(g^-2s)^tO-i)(s) + t(j+i)(s), j > 1. (4.18) 

Using the variable in (4.12), and normalizing \.^^\s) by 

rO)(t) = (,0-i)g^^-i^^)^tW(.), t = qs\ (4.19) 

one can write the fusion relation (4.18) in terms of t-polynomial operators T^^\t): 

rO) (i)r (2) {uji-H) = h{uj^-H)^h{uj'^+^-h)^r^^-^^ (t) + r^^'+i) (t) (j > i) (4.20) 

with r(°) = 0,T(i) = /, and h{t) in (2.14). 

We now discuss the algebraic Bcthc ansatz of XXZ chain of spin-^^ [44]. With the pseudo- 
vacuum J7 in (4.17), A{s)Q = a{s)^fl , D(s)r2 = a{q'^~^s)^^l , and C{s)Q = 0. For m square-distinct 

L 

non-zero complex numbers sj, 1 < j < m, we consider the vector $(si, . . . , Sm) = YliLi B(si)0 
C*^ with = M^zMzHlIi^ and define the s-function 

m m 

A{s; si,...,sm) = a(s)^ Ui + "(g'^-'s)^ /m> (4-21) 

i=l i=l 

where fg ,,' are functions in (4.6). The regular-criterion for the s-function (4.21) in the non-zero 
s-domain is given by the Bethe equation 

-$^ = -cz^™n4^^' i = l,...,m, (4.22) 

in which case by (4.7) and the same argument in Proposition 3.1, $(si, . . . , Sm) is an eigenvector 
of the transfer matrices t(,s) with the eigenvalue A(s; si, . . . , Sm)- Using the variable t in (4.19) and 
ti = qsj, (4.22) and (4.21) now take the form 

— ^ -«-na-.-*): (-3) 

r L(d-2) J r, r F(uj~^t) L(d-l)-2m , , Fiujt] 

s^gTA(,; ,1, . . . , Sm) = Hu'-H)^^^ + cj^^hiu-H)^^. (4.24) 



In the root-of-unity case where both q,uj are A^th roots of unity, we are going to derive the 

2 



evaluation polynomial for the root-of-unity symmetry study of the XXZ chain of spin-^y- with 



2 < d < N . Define the i-polynomial 



I d-2 



H{t) (= HS)) := = n (1 - ^'0, 2<d<N. (4.25) 
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Using the above function H{t), one writes the Bethe equation (4.23) in a similar form as (3.6): 

-im^'-^ \d ('"^ nt):=n(i-tr'*). (4.26) 

where 0<r<A^— l,r = -^('^ ^ ^™ 5"^) (mod A^). Now we consider only those Bethe states 
$(si, . . . , Sm) so that {ii, • • • , i^} contains no complete A^-cyclic string, (which means {w-'toljgZj^f 
for some io 7^ 0). Normalize the T^-') -operator (4.19) by^ 

^-(i- i)g-(l_ t^)TO-)(^) 

rAf+j-3 

By (4.24), the T*^^)-eigenvalue A*^^^ (t; ^i, . . . , t^) for the Bethe state $(si, . . . , Sm) is the t-polynomial 
characterized by the relation 

A(2)(t;ti,...,t^)F(t) = uj-^H{t)^F{u-^t) + H{u-H)^F{ut). (4.28) 

The relation (4.20) in turn yields the fusion relation of T^^'^'s with T(°) = 0, T^'^\t) = H{u}'^t)^ 
and 

T0)(t)T(2)(a;^-4) = Lo-^' H{LO^-H)^J^^-^\t) + H{u^-H)^J^^^^\t), j > 1. (4.29) 

Using the same argument as [42], (where formulas (3.8) and (3.9) correspond to (4.28) and (4.29) 
here), one obtains the expression of T^-'^ -eigenvalue X^^\t; ti, . . . , t„i) for the Bethe state $(si, . . . , s^): 

j — 1 Yj / k—lj.\L —kr 

which is a t-polynomial by the Bethe equation (4.26). Then one arrives the boundary fusion relation 
for T^-'^'s ( on the space spanned by all eigenspaces of the Bethe states): 

T(^+i) (t) = a;-^^T(^-i) {ut) + 2i? (a;-^)^. (4.31) 

Formulas (4.29) and (4.31) are the complete fusion relations of the XXZ spin chain of spin-^^ 
at roots of unity, which in d = 2 case, are equivalent to (4.13). By (4.30), one writes the T^^)- 
eigenvalue in the form 

AW(i;ti, ...,tm)= u-'fFiu-HfPevit''), 
where P6v(0 is the function defined by 

The Bethe relation (4.26) for tj's is the polynomial-condition for Pqsti in which case by (4.30), the 
fOO-eigenpolynomial and Pgv-polynomial are related by 

T(j)(t) + lj-^'-J^^-^XujH) = {io^Hy F{uj-H)F{u^-H)P6y{t^), < j < N, 
TW(t) = {Lu-HyF{Lu-HfPeY{t^). 



^For (1 = 2, the T(^) here differs from T^^^ in (4.12) only by the factor w'^-^^^": T^^^t) = co'^^-^^^'' T'^\t), hence 
in sectors S" = (mod N), it defines the same operator which was used in the discussion of [42]. 
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The above relations reflect the QQ-relations in the functional-relation setting of XXZ spin chain of 
higher spin for a proper Q-operator which encodes the root-of-unity property; the construction of 
such Q-opcrator is not known yet except the spin-i case [42]. Note that for d = 2,N, H{t) = 1 — t , 
PeviO is the Drinfeld polynomial in [18, 19, 42, 36] for the root-of-unity XXZ chain of 
spin-^, ^^^^ respectively. For a general d, the Pgv-polynomial is indeed the evaluation polynomial 
for the root-of-unity symmetry of XXZ spin chain of higher spin, which will be verified in the next 
subsection. 



4.2 The s/2-loop-algebra symmetry of in the root-of-unity XXZ spin chcdn of 
higher spin 

In this subsection, we arc going to show the s/2-loop-algebra symmetry of the XXZ chain of spin- 
with the A^th root-of-unity anisotropic parameter q, even chain-size L, odd A'', and the total spin 

= (mod iV) for 2 < d < N. The Pev-polynomial in (4.32) will be shown as the evaluation 
polynomial, i.e., the Drinfeld polynomial, for the sZ2-loop-algebra representation. For d = 2, the 
result is known by works in [15, 19, 18] , and for d = N, the conclusion is obtained in [36]. Here 
we are going to derive the root-of-unity symmetry of XXZ spin chain of higher spin along the line 
in [15, 19]. 

Define 

L ^ ^ L ^ 

= y^K2 . . . 0K2 0e±0K~ (g) • • • (g) K~ , = V K~ g) • • • g) K~ g)e+g)^2 (g • • • (g K 2 . 

'=1 i-1 L~i »=1 i-1 L-i 

The leading and lowest terms of entries of (4.16) are given by 



A± = lim.._(±.)^^,^A(.), B± = lim.._(±.)^(^-i),^^^-^^, 

C± = lim..^^(±.)^(--),^^^^-S^, D± = lim.._(±.)^-,^D(.). 



U. D_u = lim.+ ,„r±.s^Tifl^^^^Df.sV 

q-q 

Using (4.15), one finds 

q^^' =A^ = D±, T- = B+, S- = B_, S+ = C+, T+ = C_, (4.33) 

which give rise to the representation of Uq{sl2) on (g (by k^^ = ki = q^^"" , bq = T~,/o = T~^, 
ei = , fi = S~). Consider the normalized nth power of operators S^, as in [15, 36], S'^^") = 
^ rp±(n) _ ^^^j (n > 0) for a generic q, then take the limit as q being the A^th root of unity, 
where [n] = '^^Zq-i > W = HiLiH [0]! := 1. For q^ = 1, one obtains the operators. 



0<fc,<d, ki+- + kL=N [k{\\-[kL\\ '^«=1 
± ^ ' - 2^0<ki<d, kl+-+kL=N [ki]\-[kL]\ ^»=1 * j±' 



with the relation T^^^^ = RS^^^^R ^, where R is the spin-inverse operator of (g C^. Note that 

by g)f=iFj~^^^<*~^^>*-"''ef4 = g)f=ief4Fj~^^^<'~^^>'-"'', the order of powers of Ki and ei± 
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in the above formulas can be interchanged for the same operator. The ^ and S'^^^^T'^^^^ give 
rise to a sZ2-loop-algebra representation with their relations to Chevalley generators as follows: 

-Ho = Hi = ^- Eo = r-(^), El = Fo = T+(^), Fi = S'^^'l (4.34) 

The above Chevalley basis is related to the mode basis of the sl2-loop algebra, e{n), f{n), h{n) {n G 
Z) by e(0) = 5+(^),/(0) = e(-l) = T+W,/(l) = T'W, ^0) = ^. 

JV+l-2i 

Set m = N, Si = xq 2 [i = 1, . . . , N) in (4.7), one obtains 
A{s) nili B(.,) = nili B(5,)A(.) + ill^^Bis) UI2 B(.OA(.i), A, B ^ D, C; 

D(.) nili B(5,) = nili b(.,)d(5) - ^;gfiyff-, B(.) uli' Bis^sN), D, B - A, c. 

Multiplying each B(sj),C(si) in the above formulas by the factor (±s)j g 2 ^qr_g-J^)-J^ 
at a generic q, then taking the 00-limit for sf^ at = 1, one arrives 

A(5)5±w =5±(^)A(s) +(T5g'^)^'5±(^-l)^^gT^^ 

D(s)5±W =5±(^)D(s) -(Tsg'^)±^5±(^-l)^gg±^^ 
A(5)rT{iV) =rTWA(s) -(^sg^)±l^^(^-l)?ggT^^ 

D(s)rTW =rTWD(s) +(Tsg'^)±l^T(^-l)g^'jg±^^ 

which imply 

[t(.),5±(^)] = (.g^)=^^5^(^-i)g:;(,^^-g-^^), [t(.),r^W] = -(.g^)^^r^(^-^)^g(/^-g-^^), 

hence S^^^\T^^^^ commute with the transfer matrix t(s) in sectors = (mod N). Denote 

by V the subspace of (g) consisting all vectors with = (mod A'"). Therefore the T^^)- 
and T^-'^-operators (4.27) when restricted on V possess the s^2-loop-algebra symmetry by (4.34). 
Hence each degenerate eigenspace gives rise to a finite-dimensional irreducible representation of 

s/2-loop algebra. It is shown in [12] that every such representation is obtained by an irreducible 

M 

representation of © s/2 on <E)^Li C^*" through evaluating the loop-variable at a finite number of 
distinct non-zero complex numbers, ai, . . . ,aM, with the evaluation (Drinfeld) polynomial defined 
by P{^) = YlkLii^ — fflfcO^*'"^- Hence the representation space is generated the highest weight 
vector \v), i.e., the unique vector (up to scalars) with the highest weight among ^^-eigenspaces, 
and P(^) is expressed by 

q+{N)rrp-(N)r 

pio = EA^^(-o^ — 7-K-2 — \v) = f^r\v) 

with the degree equal to Efeli(4 - 1) ([18] (1.9), [19](1.17)). The polynomial P(^) can be deter- 
mined by the following current ([19] (1.20))^^^ 

00 

E-iO = J2fin)e, (4.35) 

n=0 

^''Hcrc wc use the current slightly different from the one used in [19] (1.20) since we employ another, but equivalent, 
set of representatives for the Chevalley basis in this paper. 



24 



whose pole-structure coincides with the zero-structure of -P(C)- Conjecturally, every Bethe state is 
annihilated by S+W,r+W ([19] (1.11), [18] (5.1) (5.2) and Sect. 5.2), by which a Bethe state 
$(si, . . . , Sra) is the highest weight vector of the s/2-loop-algebra representation to which it belongs. 
We are going to determine the current (4.35) by the "vanishing A^-string" method in [19]. 
Consider the average of a commuting family of operators 0(s): 

(0) (= (0)(.^))= n OiQ's). 

1=0 

As in (2.25) of the generalized r*-^)-model, one can determine the averages of monodromy matrix 
(4.16) of the spin-^^ six vertex model as follows: 

Proposition 4.1 For a positive integer L (no even condition required), the average of the Lth 
monodromy m,atrix (4.16) for 2 < d < N coincides with the classical Lth monodromy associated to 
(4.14); as a consequence, the averages of the entries are given by^^ 

(A) = (D) = (s^ - s-^)^, (B) = (C) = 0. (4.36) 

Proof. Parallel to (2.32) in the generalized r^^) -model, the following relations about averages of the 
(L + 1)- and L-th monodromy matrix hold: 



(Al+i) = (Al) (Lo,o) + (Bl) (Li,o) + Ek=i [^L]N-k{s)[BLUq^'^s) 'i^) 

(Af4 
-k,N 



{Bl+i) = {Al) ® (Lo,i> + (Bl) ^ (Li,i) + Ek=i'[^L]N-kis)[BLUq''-''s) lf^'\s), 



Al+1, Bl+i,Al, Bl — > Cl+i, D^+i, Cl, Dl, 

where [0]„(s) := Y[7=o ^i^^^) operators 0(s) and non-negative integers n. When L = 1, the 
average of monodromy matrix is given by (Lo,o) = (Li,i) = — s~^ , and (Lo,i) = (Li^o) = 0- Hence 
it suffices to show dj^(^^\s) = U^f^^\s) = ioi 1 < k < N — 1. Since Lo,o, are interchanged 

under the conjugation of the spin-inverse operator of C^, and the same for Lo,i, Li^O; the definition 
of and L^jf^^^^ in turn yield the relation: Lj^^"^]y(s) = R^^L^jf^Q^^\s)R. So one needs only 

to show = for I < k < N — I. Note that Uj's are d x d matrices with Lj^j diagonal, and 

^i. j 7^ i) upper- or lower-triangular with Lfj = 0. This implies l[,^^^^(s) = for k > d. Hence 
the first relation in above about the averages of A^+i and Aj;, for L = 1 yields 

(A2) = (Lo,o) ^ (Lo,o) + X:Vo,o]iV-.(^)Lo^ ^ iZ^'Hs). (4.37) 

k=l 



Since Lq for 1 < A; < cZ — 1 are linear independent lower-triangular matrices, the invariant of (A2) 
under s 1-^ qs implies that the same holds for each term in the summation of the above formula. 
Hence [lo,o]N-k{s)l^o,i ® ^IT'^i^) = [lo,o]N-k{qs)l^o,i » ^k^'^ (qs) , then 

Ln.o(.s)Lf[i ^ lE+'V.s) = lo,o{q''-'s)l'o,i ® Li7'^(«^)' ^<k<d-l. 



The (B) = for d = 2 is the formula (1.36) in [19] in the vanishing discussion for the complete iV-string Bethe 
ansatz. 
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Compare the {k, 0)th and (d — 1, d — 1 — k)tYi entries of Lg^ in the above equality, then one finds 

where a{s) is the function in (4.2). Hence ^^jy^L^^^^^ (s) = ^l^gf-i-ll^ L^^"^^\s), which by ^^^^ ^ 
a{q^~'^-Ks) ' i™pli6S l[,^^^-'(s) = for 1 < A; < c? — 2. For k = d — 1, one has the relation 
a{q'^-h)l'al^^o\s) = a(s)Li^+o^(gs). Since l^^tfi i^) = V'(s)L?,o^ where V'(s) is a (/j-term sum of 
(TV— (i+l)-productsof a((^'s)'s, one has a(g'^~-^s)V'(s) = a{s)ip{qs). Therefore, tlj{s) = ceYl'lz^a{q'^s) 
where a is a function invariant under s i— > gs, hence being a constant scalar by the degree consid- 
eration . By this, the relation (4.37) becomes 

N-1 

(A2) = (Lo,o) <8 (Lo,o) + a( n «(«'^))Lo^' ® ^ffi^^ « ^ C. 

i=0 



From the commutation relation of A and B in (4.9), (A2) commutes with 62(5) = (Lo,o "X* Lo,i + 
Lo,i® Li^o)('S). Since [62(5), LQ^^(g)L^Q^] / 0, the scalar a is equal to 0, which implies L^^^'o'*(*) ~ ^• 
Therefore L^^'''^'*(s) = for all k, hence follows the result. □ 

We now relate the sZ2-loop-algebra generators (4.34) with the vanishing averages of B, C in 
(4.36). In later discussions, we shall consider the logarithmic derivative of relations in (4.5). 
For simplicity, we shall use the subscripts of variables s,q, ... to indicate the logarithmic partial- 
derivative sds, qdq, ... of operators or functions, c.g, = s{dsB), Bg = q{dqB), {fs,s')s = s{dsfs,s'), 
(/s,s')s' = s'ids'fs,s') etc. By (4.36), the s-derivative of (B),(C) vanishes, and the g-derivative of 
(B), (C) is expressed by 

N-1 N-1 

(B), = 0, (B), = ^ B,(sg") n B(sg^), B ^ C. 

n=0 i=0,j^n 

By (4.33), the leading and lowest terms of (B)^, (C)^ are given by 



q-{N) _ lim,^o(-s)^(^-l)(B),(s) rp-(N) _ lime^ooS-^(^-i)(B)g(5) . 

~ 2iV(l-a;)-(l-a;^-i) ' ~ 2N{l-u))-{l-ui'^-'^) ' 

q+(N) _ lim,^oo^-^(-^-i)(C)<,(8) rp+(N) _ lim,^o(-s)^(-^-i)(C),(s) 

~ 2Ar(l-a;)-(l-w'^-i) ' ~ 2N(1-oj)-(1-uj'^ -^) ' 



(4.38) 



2N{l-u))-{l-ui^-^) ' 2Ar(l-a;)-(l-a)^-i) 

As in (2.1)-(2.14) of [19], for a given function (p{s), there associates the current: 

BW(s) = BS^)(s) + Bf)(s) (4.39) 



invariant under s ^ qs, where B['' \s) = (B)„ B'f'is) = En=oV(s9")B.(sg") a=o,W B(sg»), 
which satisfies the relations: [BW(s), S(s')] = [BW(s), BW(s')] = ([19] (1.39) (1.40)). For a 
Bethe state $(si, . . . ,Sjn), we are going to determine the Fabricius-McCoy current, which means 
the current (4.39) for a suitable ip{s) with the following property: 

m' m' 
t(s)(n BW(xO)$(si, ...,sm) = {l[ BW(xO)t(s)$(si, ...,Sm) (4.40) 
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for all xi and integer m' >1. By differentiating relations in (4.7), one has 

Hs)Bs{Sn) mi,^n B(S.) = (n™ 1 fs,)Bs{SnmT=l,^n B(5^))A(s) 

- Er=i {9s,k UT=i,^k fk,)^^ B{smT=i,^k B(5.))A(sfc) + iUT=i fs,i)s^ {I{T=i ^{si))k{s) 

(A(.) mi B(s.)), = {uT=i fs,) (nrii b(..)a(.)), - etH {9s,k uT=i,^k a,)^ B(.)(n^i,^, b(.,))a(.,) 

+ (n^i fs,i\ (IlT=i B(s.))A(.) - ET=i 9sAUT=i,i^k h,i) (b(s) Y[T=i,i^k ksmsk))^ ; 

^{s), fs,ij fk,ij 9s,k ^{s), fi,s, fi,ki 9k,s- 

Set m = N, Si = xq^ for 1 < z < N{= 0) G Zjv in the above relations, then using the relations 

fk,k+l = 0) ni=0 fs,iYli=0 fi,s = Ij {I\i=o,ij^k fk,i)q = (ni=o,i^fc = q-q-^ ' 
{Y\i=0^i^k fk,i)s„ = (<^fc,n ~ Sk,n~l) q-q-^ i (ni=o,i^fe /i,fe)sn ~ {^k,n ~ ^k,n+l) q_q-L , 

one finds 

/K{s)Bf\x) = BfHx)A{s) - -^E^roi5,,„B(.)(niI-o,UB(xgO)A(xg"), 
D{s)Bf\x) = Bf\x)D{s) + ^E:!=o' 9s,nB{sml-^^^^^^ 

A(.)Bf = Bf )(x)A(s) + _|^E^=o'5.,nB(.)(9.(zg"+i) - v^(xg"))(niIo,W B(xg»))A(xg«), 

D(s)Bf)(x) = Bf)(x)D(s) + ^En5s,nB(s)((^(xc?"-i)-v^(xg"))(niIo,UB(^^^ 
Hence 

A(s)BW(x) = BW(x)A(s) + |^E^roi5,,„(c^(xg«+i) - ^(zg") - l)(niIo,UB(xgO)A(xg-); 
D(.)BW(x) = BW(x)D(.) + ^ EtV ^..nl^'N'^-^) - + l)(n£o,W B(xgO)DN"). 

The relation (B) = yields 

(a=o,U B(xgO)A(xg")cl.(si, ...,sm) = (OilaU B(xgO)(a^(xg") D^i /.g",.J^(si, • • • , ^m), 
(niloj^n B(a;gO)D(xg")$(si, . . . , s„) = (n.=ol^n B{xmaHxq-+''-^) UZi . . . , s^), 

hence the condition (4.40) for m' = 1 is provided by the following constraint of (p{s): 

m m 

(ifisq) - ip{s) - l)a^(s) n + {^{sq-') - v{s) + l)a^(s/-i) fs„.qn = (4.41) 

1=1 i=i 

for all s. Use the variable t in (4.19) and write the function (4.32) by Pev{t^) = Ylk=o Pi^^^) 
with p{t) := pf^pl'l) ■ Up to the factor q-^'^-'^^^+"'t^+''F{ujt)F{LO-H)U.tZo{^ - ^"t)'^, (4.41) is 
equivalent to (99(sg) — (/3(s) — l)p(a;^^t) = {ip{s)—ip{sq~^) — l)p{t). Hence (p{sq)—ip{s) — l = a{t'^)p{t) 
for some function a{t^). The relation (p{sq^) = (f{s) in turn yields a(t^) = jv — ~) hence the 
(^-condition (4.41) is equivalent to the following equation of ip: 

^W-^W-l = ^. (4.42) 

Up to additive s^-functions, the above equation has the unique solution given by 

o,(.) = ^^&°M^ .(t) = ^^^ (4 43) 

^^'> Pevit^) ' ^^'^ F{t)F{ujty ^ ''^ 
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Note that by {B)s = 0, all solutions ip{s) of (4.42) define the same current (4.39). Therefore 
we have shown that with (p in (4.43), the relation (4.40) holds for m' = 1, then also true for all 
positive integers m' by the same induction argument in Sect. 2 D of [19]. Hence we obtain the 
Fabricius-McCoy current given by (4.39) with (fi{s) in (4.43). Now we show the following result. 

Theorem 4.1 Let $(si, . . . , Sm) be the Bethe state for the solution, ti (= Qs|), 1 < i < m, of the 
Bethe equation (4.26), and Pev{Cj he the polynomial (4.32) defined by ti's. Then the current -E~(0 
in (4.35) and the Fabricius-McCoy current B^^\s) in (4.40) are related by 

(_s)M-f'-i)BW(s) = 2iV('n (1 - u'))E-{t''). 

i=0 

As a consequence, j^Peyi^) is the evaluation polynomial for the sl2-loop-algebra representation 
determined by ^{si, . . . , Sm)- 

Proof When acting on the Bethe state $(si, ...,Sm), the i^-operators, E-{t^) and (-s)^(^~^)b(^)(s), 
both satisfy the property (4.40), hence they are proportional to each other. By {B)s = 0, the low- 
est term of {—s)^^^~^^B'^^\s) is contributed only by the term from B^''^^(s), which is equal to 
2Arn^^-i(l - u;^)^-^^) by (4.38); while the current E-{C) has the constant term E-{Q) = /(O) = 
5'-{^), Hence follow the results. □ 

Remark. By the same argument in Sect. HI A of [19], the current B*^^) (s) in the above proposition 
has poles only at the zeros of PQY{t^), which is consistent with that for E~{^). □ 



4.3 Comparison of the XXZ chain of spin-^^^ at Nth. root-of-unity and super- 
integrable CPM 

It was noticed [1, 5] that the superintegrable CPM and the XXZ chain of spin-^^^^ for = 1 share 
the same Bethe equation (up to the phrase factors). In this subsection, we compare the symmetry 
of the degenerate eigenstates between two models, and provide the answer to a question raised in 
[36]. 

Under the change of variables t = utt' , the polynomial (4.25) for the XXZ chain of spin-^^^^ 
and (3.11) for superintegrable CPM are related by 

HN{t) = Hcp{t'), 

by which, a six-vertex-Bethe-solution {tj}™ of (4.26) is equivalent to a CPM-Bethe-solution {t'^^^ 
of (3.10) for sectors, 5|v = (Pa + -P6)cp = r (mod iV), the T -eigenvalue X'^'^\t;ti, ...,tm) in 
(4.28) the same as T^'^\t') in (3.13), and the identical evaluation polynomials, (4.32) and (3.14) : 
Pgv(i'^) = u}~^'^^''Vcp(t'^ ) ([36] Proposition 4.1). In the theory of Onsager algebra symmetry of 
superintegrable CPM, the polynomial VcpiO is well-understood as follows. It is known that the 
Onsager algebra can be realized as the Lie-subalgebra of the sZ2-loop algebra fixed by a standard 
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sZ2-involution and inverting the loop- variable [38]. The (finite-dimensional) Onsager-algebra rep- 
resentation was known in [14] , and the rigorous mathematical theory has been fully developed and 
understood in [13], in particular, the Hermitian irreducible representations have been completely 
classified. By this, one finds VcpiO ^ simple ^-polynomial with only negative real roots and the 
degree tue = the integral part [ ^^ ^ ^"' J, which implies the dimension 2™^ for the degenerate 
r^^^-eigenspace as an irreducible Onsager-algebra representation ([40] Theorem 2 and formula (70)). 
Correspondingly in the XXZ spin chain, one concludes only the spin-^ representations occur in the 
si2-loop-algebra representation generated by the Bethe state $(s, si, . . . , Sm) with 2™^ -degeneracy 
for the t'^^^-eigenspace, (by which, follows the answer of the question raised in [36] Proposition 
4.2). Prom the above discussion, the root-of- unity XXZ chain of spin-^^^^ and superintegrable 
CPM share the equivalent Bethe equation and evaluation polynomial, with the same degeneracy 
of eigenstates for certain sectors, but they carry the different type of symmetry structure. The 
Onsager-algebra symmetry of the superintegrable CPM is inherited from the Baxter's Q-matrix, 
indeed generated by the symmetry operators of the quantum Hamiltonian chain (see [40] Sect. 
3); while the root-of-unity symmetry of the XXZ spin chain arises from the g-derivative (4.38) of 
vanishing-avcragc-cntrics of the monodromy matrix. It is pertinent to ask whether there exists a 
larger symmetry algebra than the Onsager algebra in the superintegrable CPM as suggested by the 
XXZ spin chain. But, the answer has not been found yet. The algebraic-Bethe-ansatz discussion of 
r(^^-models in section 3 could possibly provide certain clues to this end, though not clear at present. 
This is also one of the reasons that we conduct the algebraic-Bethe-ansatz study of r^^^ -model in 
this work. 



5 Concluding Remarks 

In this article, we have made a systematic account on the algebraic approach about the fusion op- 
erators of the generalized r'^^^-model using ABCD-algebra techniques. The recursive fusion relation 
(2.23) follows automatically from the construction of fused L-operators. Under the modest and 
seemingly innocuous conjecture (2.34) supported by computational evidence in cases, we produce 
a logical explanation about the validity of boundary fusion relation (2.31) for the generalized r^^^- 
model, by which one can use the separation-of-variables method to solve the eigenvector problem 
of the model in generic cases ( [24] Theorem 2). On two special classes of r^^^ -model centered 
at the superintegrable one, we perform the algebraic-Bethe-ansatz technique to study the Bethe 
equation and Bethe states. The efforts enable us to reconstruct the r^^^ -eigenvalues of certain 
sectors among the complete spectra previously known in the theory of superintegrable CPM [1, 5] 
through the functional-relation method. With the similar argument, we discuss the root-of-unity 
XXZ spin chain of higher spin, and obtain the fusion relations with truncation identity as in the 
spin-i case. Furthermore, the algcbraic-Bcthe-ansatz technique can produce better quantitative 
results in the root-of-unity XXZ spin chain about solutions of the Bethe equation than those in the 
7-(2)_j]2odel partly due to the symmetric form of the six-vertex i?-matrix. Indeed from the Bethe 
relation and fusion relations, we successively extract the evaluation polynomial for the sZ2-loop- 
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algebra symmetry in the root-of-unity XXZ spin chain of higher spin, verified in subsection 4.2 
following the line of the spin-^ case in [15, 19]. Since results in the root-of-unity symmetry of XXZ 
spin chain bear a remarkable quantitative and scmi-qualitativc resemblance to the Onsager-algebra 
symmetry of superintegrable CPM, it would be desirable to have a Baxter's Q-operator for the 
XXZ spin chain of higher spin at roots of unity to pursuit the functional-relation study about the 
symmetry as in the spin-i case [42]. Such Q-operator has not been found yet. The aid of a special 
Q-operator in the XXZ spin chain to encode the root-of-unity symmetry will provide a useful tool to 
make a detailed investigation about the degeneracy, much in the same way as the Onsager-algebra 
symmetry does in the superintegrable CPM. In particular, in the spin-^^^^ case it is possible to 
construct the Baxter's Q-operator of the root-of-unity XXZ spin chain so that the comparison in 
subsection 4.3 will be understood in a more satisfactory manner. A program along this line is under 
consideration, and progress would be expected. In this work, we discuss the symmetry about the 
generalized r^^' -model, and the XXZ spin chain of higher spin. We also hope that our results will 
eventually lead to the understanding of other models, such as the root-of-unity eight-vertex model 
in [16, 17, 21, 22]. This programme is now under progress. For the root-of-unity symmetry of 
XXZ spin chain , our discussion can be applied to a more general setting. But, just to keep things 
simple, we restrict our attention in this paper only to even L, odd N and 5^ = (mod A^), and 
leave possible generalizations, applications or implications to future work. 
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